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We discuss the fluctuation properties of diagonal matrix elements in the semiclassical 
limit in chaotic systems. For extended observables, covering a phase space area of many times 
Planck's constant, both classical and quantal distributions are Gaussian. If the observable is 
a projection onto a single state or an incoherent projection onto several states classical and 
quantal distribution differ, but the mean and the variance are still obtainable from classical 
considerations. 

§1. Introduction 

A central tool in the semiclassical analysis of chaotic systems is the Gutzwiller 

n n 

trace formula which relates quantal properties to classical periodic orbitsEP'EP. The 
extensions of Gutzwiller's trace formula that include matrix elements in the density 
of statescP'o) allow to relate distributions of quantum matrix elements and distri- 
butions of classical trajectory segments Ef' 13. Specifically, quantum fluctuations are 
related to fluctuations in ensembles of classical trajectory segments of length Heisen- 
berg time. As the classical limit is approached, the Heisenberg time diverges, and the 
fluctuations vanish. This approach can therefore be used™ to study the approach 
to ergodicity in quantum wave functions in chaotic systems. Integrable systems are 
not ergodic on the energy shell and therefore the fluctuations do not vanish in the 
semiclassical limit, as shown by MehligQ). The form factors for the matrix element 
weighted densities of states have an additional constant contribution tf and different 
behaviour in their rigiditvtP. Moreover, these relations can be used to analyse the 
effects near bifurcations B and to study localisation effects in wave functions. Be- 
yond the static properties these can also be used to analyse the behaviour of cross 
sections©©. 

The aim here is to analyse further an observation first made in the calculation 
of cross sections. By Fermi's golden rule photodissociation cross sections contain 
matrix elements for observables that are projections onto initial states. In this 
limit semiclassical arguments become suspect since the smoothness of the observable 
that was assumed in the derivation of the trace formulas is no longer guaranteed. 
However, numerical experiments indicate that even in this limit at least the widths 
of the distributions can be obtained from classical considerations. This behaviour of 
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matrix elements of projection operators will be investigated further. 

The behaviour of matrix elements was studied in a variety of systems, including 
hydrogen in magnetic fields B' S 'E3), the bakers mapS 1 , a perturbed cat mapli^, and 
various billiardsL7'0). We here use the kicked standard map as the main example llH. 
The main advantage of maps is that they allow to study the semiclassical regime more 
easily and more deeply than continuous systems. 

The outline of the paper is as follows. In section ^ we summarise information 
on the quantised standard map, on the semiclassical arguments that relate classical 
and quantal matrix elements and on the behaviour expected from random matrix 
theory. The numerical and analytical results on projectors are presented in section 
||. We conclude with a summary and an outlook in section 

§2. Preliminaries 

2.1. The matrix element weighted density of states 

In& we derived formulae for the width of the quantum matrix element distribu- 
tion in smooth systems and proposed a formula for maps by analogy to the smooth 
results. To set the stage for the present investigation and to explain the origin of 
this formula, we here sketch the equivalent derivation for quantised maps on a finite 
phase space (similar to Carvalho et alE3). 

Let U be the unitary operator for the time evolution, A the observable and 
Nd the dimension of the Hilbert space. Then there are Njj eigenphases <j)a and 
eigenstates so that 

0- = X>>^ M </*l- (2-1) 

The matrix elements = (/j,\A\fi) are conveniently contained in the weighted den- 
sity 

p A (0) = £A^>-^)- ( 2 - 2 ) 
it=i 

It will be useful to work not with pure delta functions but with smeared out objects 
so that products are well defined. Rather than the Lorentzians or Gaussians used 
in previous works, the relation to finite trigonometric sums suggests to work with 
sine- functions, i.e. 

_ J_sin[(2/V + lW21 _ l_ £ 
v ' 2tt sm(x 2) 2vr ^ M V ' 

v ' / n=—N 

of width e ~ 2ir/(2N + 1). The advantage of this function is that the smoothed 
density can be expressed as a finite sum of traces of powers of the propagator, viz. 

1 N Nd 1 N 

Pa{^) = — £ £V^-*) = — Y, tr(AU n )e-^. (2-4) 

n=—N ju=l n=-N 
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If N is significantly smaller than Njj some of the sinc-functions will overlap and not 
all individual eigenphases can be identified. Very large ./V are impracticable from the 
point of view of a semiclassical approximation to be made below. Thus a natural 
choice is N of the order of Njj so that all eigenstates can be resolved in the mean, 
yet the semiclassical approximations are still reasonable. 

The autocorrelation function of the smoothed density of states becomes 

Ca(0) = ±- d<t> PA {<t> - 9/2)p A (0 + 9/2) 
2tt Jo 

4vr2 



1 N 

— 2 £ |tr(AC/")|V" e (2-5) 

n=-N 



Up to some factors, the Fourier transform of this correlation function is the form 
factor. Here, since the spectrum is confined to the unit circle and a periodically 

continued interval in </>, the form factor depends on a discrete argument only. We 

(A) 

split off the density of states and define the form factor Kn by 

= — £ K™**. (2-6) 

n=— oo 

When calculated with the smoothed density of states only the first few elements up 
to Nd are included, but they do not differ from the exact ones, 

K^ = ^-\tv(AU n )\ 2 ■ (2-7) 

As in & , all the information needed on matrix elements is contained in the semiclas- 
sical evaluation of K n . 

Within the semiclassical approximation the trace is calculated from a station- 
ary phase approximation to the propagator, which leads to an expression involving 
classical periodic orbits. Specifically, 



tv(AU n )\ sc = J2 A 

7G(n) 



where the sum extends over all distinct orbits 7 which close after n iterations. S~ is 
the action associated with the trajectory, w-y the amplitude including phases and Ay 
the sum of the observable A over all points of the trajectory. Of the amplitudes we 
need that their modulus squared equals the classical periodic orbit weight, |u> 7 | 2 = 
1/1 det(l — M 7 )|, where M 7 is the linearisation perpendicular to the orbit. Thus the 
semiclassical form factor becomes 

\tv (AU n )\ sc \ 2 = £ A 1 A*,w 1 w*,e i{s ->- s t' /n g ^ |A 7 | 2 |w 7 | 2 . (2-9) 

7,7'G(n) 76(n) 

The neglect of the contributions from orbit pairs (7,7') with 7/7' constitutes the 
diagonal approximationllB'. The factor g arises from the possibility of a quantum 
interference of classically distinct orbits with the same action, e.g. in the case of time 
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reversal invariant systems, where a path and its time reversed image correspond to 
different objects in phase space, but with the same action. Because of the quantum 
interference between amplitudes, it contributes with twice the classical weight, hence 
g = 2. In systems without such symmetries, like the unitary ensemble, g = 1. 

The extraction of the variance of the matrix elements from the form factors 
proceeds as in Ref. 0). The contributions Ay from orbits with similar period consist 
of two parts, one proportional to the mean of the observable and to the length of 
the orbit and one capturing the fluctuations, i.e. (A 2 ) = (Aj) 2 n 2 + an. Taking the 
square, averaging over orbits of similar length and exploiting the periodic orbit sum 
rule (which absorbs one factor of n) gives 

K - A) ~ 9{A)2 w D +9a ik (2 ' 10) 

The linear increase of the first term is the small n expansion of the form factor, correct 
up to n = Nd for the unitary ensemble but approximate only for the orthogonal 
ensemble© . If the mean of the observable vanishes is already the variance 



of the matrix elements and thus contained in the second term of ( 2-10 ). Since the 
periodic orbit sum rule is an asymptotic statement, we take n as large as possible, 
but limited by the requirement that the diagonal approximation be reasonable and 
a semiclassical approximation possible. We therefore take n = Njj, the Heisenberg 
time. 

This calculation relates only the first and second moments of classical and quan- 
tal distributions. The main aim here will be to compare the full distributions, espe- 
cially in the limit of singular observables. 

2.2. Connection to random matrix theory 

As demonstrated in a number of studies 0' III)'©, the short range statistical 
properties of quantised chaotic maps are in good agreement with those of ensembles 
of random matrices. So in the present context one can ask for the prediction of ran- 
dom matrix theory for the diagonal matrix elements. In particular, for the kind of 
maps and observables studied here, the relevant input is the statistics of eigenvector 
components. For the unitary ensembles, the distribution of real and imaginary parts 
can be calculated from the assumption of independent elements and the normalisa- 
tion condition. For sufficiently large dimensions, this results in essentially Gaussian 
distributions, 




P(x,y) = ^ e -^*%/-^r~ W . (2-11) 

If An = (/j,\A\/jl) denotes the matrix element of A with respect to the basis set, the 
matrix element with respect to an eigenstate /i described by the amplitudes x^ + iy^ 
becomes 



A = + (2-12) 
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The probability density for the distribution of the matrix elements then becomes 

P(A) = / 6(A-J2M4 + yl» II P(x^,y^dx^. (2-13) 
J p- p=\ 

By the usual manipulations this becomes a Gaussian in the limit of large Nd, 



P{A ) = jJ[E_ e ~N D (A-(A)f/2«l (2 . 14) 
27TCT 



where 

are the average and variance of the matrix elements of A in the original basis. 
2.3. The map 

We close this section with a description of the quantum version of the kicked 
rotator used here. The kicking potential 

V((f>) = fc(cos^-7sin(20)) (2-16) 

gives rise to the classical map 

Pn+i =Pn + &(sin + 27 cos(20)) 

4>n+l = 4>n+ Pn+l ■ 

The parameter k is the kick strength and 7 is used to break a reflection symmetry 
of the map. It is obviously periodic of period 2ir in <p. As we will study the quantum 
map at the case of resonance, the momentum p is not unbounded but also restricted 
to be periodic with period 2nr, where r is the number of classical resonances. The 
parameter 7 can be used to change continuously between the cases with time reversal 
symmetry (7 = 0) and without (we take 7 = 1, for intermediate behaviour seeElrlHJ). 

In the earlier quantum models of this map, phase space was taken to be an in- 
finitely extended cylinder where arbitrarily large momenta were allowed, the quan- 
tum map was infinite. However, it was noted that for rational values of Ti, a consistent 
quantisation with a finite map can be achieved. If the classical phase space extends 
from — irr to +irr in momentum and if the quantum map contains Ne> states, the 
effective value of Planck's constant has a width of 27T in angle, the total phase space 
area is 47r 2 r. This has to be divided into Nd quantum cells, hence h = Air 2 r/N or 
h = 2irr /Nd- We work in a momentum basis with wave functions {<p\n) = e mi ^ ' j\phx 
and momentum eigenvalues p n = fin = lixm /Nd ■ To avoid common factors between 
r and Nd, we take the dimension odd, Nd = 2iVi + l. The quantum unitary operator 
then has the form 

U nm = — e -f (n 2 +m 2 ) g-iV^/h-iin-m)^ (2-17) 
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where the sum extends over the points <pj = 2ir j/No- In order to have full symmetry 
over the finite momentum interval, 

one should take the number of resonances to be a multiple of 4 (if r is just even 
translation by Np produces a sign change). 

§3. Quantum and classical matrix element distributions 

In this section we compare the distributions of matrix elements in several cases. 
We consider the quantum kicked rotator (2-16, 2-17) and its classical counterpart 



(2-17). The classical kicking strength is taken in be K = kfi = 7 in order to ensure 
complete ergodicity over the whole phase spaceclP. The matrix dimension Nd = 201 
turns out to be sufficient to provide appropriate statistics. 

3.1. Chaotic delocalised states 

The observable we consider is a projection onto a subset, I(m), of the unper- 
turbed momentum basis states of length m, 

A= £ \n)(n\. (3-1) 

n£l(m) 



Its matrix elements in the basis of the eigenstates, fi, of the evolution operator fl2-17| ) 
are 

A, = ^\A\^)= £ \{n\n)\ 2 . (3-2) 

n£l(m) 

The projection covers a fraction p = m/No of the phase space both quantum me- 
chanically and classically. 

Due to normalisation of the eigenstates the mean value of A^ is {A) = p = 
m/Nij. The second moment, (A 2 ), is defined as 

^ 2 ) = ]v-E E Mk)\ 2 Ml)\ 2 . (3-3) 
D M k,lel(m) 

which for m No can be calculated using the assumptions of random matrix theory 
to give 

(A 2 )=p 2 + g^p, (3-4) 

where g = 1 (g = 2) for the case of unitary (orthogonal) symmetry. Hence the 
variance is simply a\ = gp/N . 

If m is sufficiently broad, m > 1, we may easily compare the distribution of 
these quantum matrix elements to that of the classical distribution which, in the 
limit N > m 3> 1 is essentially a Gaussian ( [2-14 , 2 • 1 5| ) . For the classical observable 



the mean and the fluctuating part can be calculated using a random-walk analogue 
due to the rapid decay of correlations between steps, 

(A) = p = m/N D , a=p(l-p). (3-5) 
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The classical fluctuations in the limit of p <C 1 are related to the quantum variance!^ 



a u N 



(3-6) 



Fig. 1 clearly demonstrates that the classical and the quantum distributions are 



orthogonal y=0 
unitary 7=1 
classical 




Fig. 1. Distribution of the matrix elements with t — NdA^. The orthogonal symmetry (7 = 0) 



has bee n rescaled according to (3-7). The continuous curve is the binomial distribution as given 
in (3-8). The parameters are: Nd — 201 and m = 20. 



indeed very similar when the quantum kicked rotator is used in the unitary symmetry 
(7 = 1). In the figure we also plotted the distribution of matrix elements for 7 = 0, 
i.e. for the orthogonal case. In the latter case we adopted a rescaling according to 
the symmetry factor ( |3-6| ) 

oi(7 = 0) = Qo\ici 



0) 

where the ^-factor (g = 2), as in ( |2-9| , pM0| , |3 4| , |3-6| ) accounts for the fact that 
in the orthogonal symmetry the orbitals and their time reversed counterparts are 
indistinguishable. Results on the transition betweenthe two ensembles can be found 
in Blumel and SmilanskyE^P and in Carvalho et al. c3 

In fact one can look at the variance of the observable and compare it to the 
mean value for different values of m/N and compare with the classical predictions. 
In Fig. 2 it is demonstrated that indeed for a wide range of m for a fixed value of 

= 201 the behaviour of the variance is the same as presented also in Fig. 1. 

3.2. Singular observables 

As the width of the projector, m, becomes smaller the assumption of an ob- 
servable that is smooth on the scale of Planck's constant is not satisfied anymore. 
Nevertheless, the data of Fig. 2 show that the variances are still in agreement with 
the semiclassical predictions. The full distribution, however, cannot remain Gaus- 
sian for m approaching one, neither classical nor quantal. On the quantum side, the 
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Fig. 2. Quantum versus semiclassical variance of A M as a function of its mean value (A) = m/No- 
The classical curve is aN (3-5). We have used Nd = 201. 



distribution becomes exponential for the unitary ensemble and a Porter Thomas dis- 
tribution P(t) = e _t//2 /\/2t for the orthogonal ensemble. The classical distribution, 
in the case of strong chaos, is a binomial distribution, i.e. the probability to hit the 
support of the observable /c-times in Nd trials is 



W k (N 



D) 



P) 



N D -k 



(3-8) 



where p = wi/Nd- From fl3-8| ) (A) = Np = m and a A = Na = Np{\ — p). The limit 
of m approaching 1 and diverging Nd is the Poisson limit of the binomial distribution: 
for p — > while Nd — ► oo so that pN = m is fixed the above distribution reduces to 



Pi- 



rn 



(3-9) 



with a clearly non-exponential dependence on k. 

Fig. 3 shows the change of the distribution of matrix elements, A^, for different 
values of m and fixed Nd = 201. The histograms nicely follow the classical distri- 
bution for large m. To highlight the differences for smaller values we show in Fig. 4 
the cases m = 1 and m = 5 on a semi- logarithmic plot. The exponential distribution 
of the quantum matrix elements is clearly visible as is the strongly non-Gaussian 
behaviour of the classical distribution. The enhanced return probability in the or- 
thogonal case and the singular behaviour of the Porter Thomas distribution give 
result in more pronounced differences between classical and quantal distributions 
than in the unitary ensemble. For instance, the distributions for m = 5 in Fig. 5 for 
orthogonal symmetry show larger differences than the corresponding ones in Fig. 3 
for the unitary ensemble. 
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Fig. 3. Distribution of the matrix elements Ay, with t = N^A^ for the case of unitary symmetry 
(7 = 1) with m = 1, 5, 10, and 20. The continuous curves are the classical binomial distributions 
We have used N D = 201. 




Fig. 4. Distribution of the matrix elements A^ with t = NdA^, for the case of un itary symmetry 
(7 = 1) with m = 1, 5, compared to the classical binomial distributions (3-8). We have used 
N D = 201. 



3.3. Matrix element distributions near bifurcations 

The previous sections have shown that while it is not always possible to obtain 
the full distribution of matrix elements from the classical distributions, at least the 
mean and the variance come out reliably. This opens the possibility to study also 
the distributions in situations where the assumptions of random matrix theory are 
not satisfied anymore, specifically in systems with mixed phase space c3 or near 
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Fig. 5. Distribution of the matrix elements with t = NjjA^ for the case of orthogonal symmetry 
( 7 — 0) with m — 1, 5, 10, and 20. The continuous curves ar e the classical binomial distributions 
(p-q). The numerical data have been rescaled according to (3-7). We have used Nd = 201. 



bifurcations of classical orbits ErEf . To illustrate the kind of behaviour that has 
to be expected near a bifurcation we consider an observable that is localised in 
a strip of width 2ir/N around the the position of the bifurcation, p = 2ir. Near 
the bifurcation the probability to stay in this box is enhanced and so will be the 
contribution of this state to the matrix elements. It is thus useful to compare directly 
the classical and quantal probabilities to return to this box. Therefore, we compare 
P c i, the classical probability for trajectories that start in the strip and return after 
./V steps and P qm = tr(AU N ) (Ref.ES 1 ). As can be seen from Fig. 6 both quantities 
show a qualitatively similar behaviour, with a rapid increase near to and a slower 
decline further above the bifurcation point. However, the maxima are shifted in a 
characteristic fashion and the scaling with N is different c3 . 



§4. Conclusions 

The semiclassical theory used in Ref. & to relate matrix elements and classical 
trajectory segments requires the observable to be smooth on scales of Planck's con- 
stant. Nevertheless, the examples discussed here show that this relation seems to 
hold all the way down to singular observables like projection operators onto single 
states. Quantum effects are always present through the symmetry factor g, but they 
are enhanced for singular observables where the form of the distribution and thus in 
particular the higher moments are different from what could be expected classically. 
Nevertheless, the possibility to analyse the second moment semiclassically should 
be helpful in analysing quantum chaos in systems with mixed phase spaces or near 
bifurcations. 
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Fig. 6. Classical and quantum return probability to the neighbourhood of bifurcation above K = 4tt 
for N D = 101 
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